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PACS 71 . 10. Li " Excited states and pairing interactions in model systems 
PACS 71 .35.-y - Excitons and related phenomena 

PACS 71 .30.+h - Metal- insulator transitions and other electronic transitions 

Abstract. - Using the projector-based renormalization method we investigate the formation of 
the excitonic insulator phase in the two-dimensional (2D) spinless Falicov-Kimball model with 
dispersive / electrons and address the existence of excitonic bound states at high temperatures 
on the semiconductor side of the semimetal-semiconductor transition. To this end we calculate 
the imaginary part of the dynamical electron-hole pair susceptibility and analyze the wave-vector 
and energy dependence of excitonic resonances emerging in the band gap. We thereby confirm the 
existence of the exciton insulator and its exciton environment within a generic two-band lattice 
model with local Coulomb attraction. 



g ] Introduction. — In solids, the Coulomb interac- 
^tion binds conduction band electrons and valence band 
'— Hioles to excitons. Normally, excitonic quasiparticles do 
^^not form the ground state but electron- hole excitations 
^that tend to decay on a very short time scale. At a 
f~ — semimetal-semiconductor transition, however, the conven- 
ed tional ground state of the crystal may become unstable 
O^Owith respect to a spontaneous formation of excitons, pro- 
^^yided the overlap or band gap between the valence and 
oonduction bands is small. Then, for low enough temper- 
^^atures, these composite bosonic quasiparticles will con- 
I dense into a macroscopic phase-coherent quantum state, 
• • thereby transforming the semimetallic or semiconducting 
. ^ configuration into an insulating one (cf. also fig. [IJ . This 
^^sb-called excitonic insulator (EI) state was theoretically 
$—1 proposed more than four decades ago [Tl[2] ; for recent re- 
. ^ yiews see [3j|4] . The EI phase realized below the critical 
temperature Tei can be perceived either as BCS conden- 
sate in the semimetal region or as Bose Einstein conden- 
sate in the semiconductor region [S]. 

The EI state is extremely rare in nature, so far there 
is no free of doubt realization in any material. At 
present, the most promising candidates are the quasi- 
2D transition-metal dichalgogenide ir-TiSe2 and the 
pressure-sensitive mixed-valence rare-earth chalgogenide 
TmSeo.45Teo.55. In lT-TiSe2, the excitonic condensate 
exerts a force on the lattice generating periodic ionic dis- 



placements [6]. For TmSeo. 45X60. 55, Hall effect, thermal 
diffusivity and heat conductivity data give strong support 
for a Bose condensed state in the pressure range between 
5 and 11 kbar below 20 K [7]. The transport anomalies 
observed at higher temperatures, in particular the strange 
increase of the electrical resistivity in a narrow pressure 
range around 8 kbar might be attributed to a free-bound 
state scattering in an exciton- rich "halo" of an EI [S] . As 
a basic prerequisite for the validity of this scenario, the 
existence of free excitons above the EI phase has to be 
proven, at least for the semiconducting region. 

The aim of this paper is to address this issue. The inves- 
tigation of Falicov-Kimball-type models seems to be min- 
imal in this respect. The original Falicov-Kimball model 
describes localized / electrons and itinerant c electrons in- 
teracting by an on-site Coulomb interaction [5]. For our 
problem, we have to allow for a possible coherence be- 
tween conduction band electrons and valence band holes 
however. This can be achieved either by including an ex- 
plicit c-f hybridization [S] or a finite / bandwidth [TU]. In- 
deed, using constrained path Monte Carlo [TT] and mean- 
field [12] techniques, the 2D Falicov-Kimball model with 
direct /-/ particle hopping has been shown to exhibit an 
excitonic ground state for intermediate Coulomb couplings 
provided that the center of the c and / bands, e'^ and , 
energetically differ. Note that around the symmetric case 
e*^ = e-^ a charge-density-wave phase is energetically more 
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Fig. 1: (Color online) Schematic finite-temperature phase di- 
agram of the EFKM as obtained from RPA [T5] and slave- 
boson [13 approaches for the particle-hole asymmetric case 
£^ ^ e'^ ^ Q < \ < . At the semimetal-semiconductor tran- 
sition the ground state of the system may become unstable 
with respect to the spontaneous formation and condensation 
of excitons [2]- The strength of the Coulomb interaction de- 
termines on which side of the metal-insulator transition the 
system is. An interesting issue is the possible existence of pre- 
formed zero-momentum excitons in the semiconducting region 
for Tei < T < Tmi, where Tmi denotes the critical temper- 
ature for the high-temperature metal-insulator transition. At 
T^ij the "Hartree" gap opens. 



Stable [TMT^. Recent Hartrcc-Fock [12 , RPA [H] and 
slave-boson |15j studies confirm this finding also for the 
3D case. In this paper we use the projector-based renor- 
malization method (PRM) [16l[T7] to calculate directly 
the excitonic pair susceptibility (up to second order in 
U). Analyzing the non-trivial frequency- and momentum- 
dependence of x{(\, uj) we are able to address the problem 
of exciton formation and condensation. 

Theoretical approach. — In order to model the 
generic situation of semiconductors or semimetals with 
short-ranged attractive Coulomb interaction between con- 
duction band (c) electrons and valence band (/) holes 
we consider an extended version of the Falicov-Kimball 
Hamiltonian (EFKM), 



number is strictly conserved jl3j . In what follows, we 
study the half-filled band case, with total electron den- 
sity (n) = {nf) + (71/) = 1. Moreover we consider a direct 
band gap situation with the maximum (minimum) of the c 
(/) band dispersion located at (tt, tt), i.e. < 0. Without 
loss of generality the c electrons are considered to be 'light' 
while the / electrons are 'heavy', i.e. \tf\ < 1, where the c 
electron hopping integral is taken to be the unit of energy, 
t"^ = 1, and e= = 0. 

The projector-based renormalization approach starts 
from the decomposition of the many-particle Hamilto- 
nian ([1} into an "unperturbed" part Hq (c and / electron 
band terms) and into a 'perturbation' Hi (Coulomb in- 
teraction term), where the unperturbed part Ti.Q clearly is 
solvable. Then, in general, TLi accounts for all transitions 
between the eigenstates of T-Lq with nonzero transition en- 
ergies. Using a series of unitary transformations to inte- 
grate out the perturbation Hi (for details see Ref. [TO] ) 
one arrives at a final Hamiltonian which is diagonal or at 
least quasi-diagonal. To evaluate the expectation value 
(A) of any operator A also the operator has to be trans- 
formed by the same unitary transformation. One of the 
main advantages of the method is to find broken symmetry 
solutions of phase transitions [12] . Note that for practical 
applications the unitary transformations should best be 
done in small steps in energy. Therefore, the evaluation of 
the transformation in each small step can be restricted to 
low orders in Hi. This procedure usually limits the valid- 
ity of the renormalization approach to parameters values 
of Hi which are of the order of those of Hq. In the present 
case, good agreement with exact results is expected for U 
values smaller than f-^ or — (cf. eq. (|3]) below). 

Excitonic insulator phase. — In a first step, let us 
address the formation of the long-range ordered EI state in 
the 2D EFKM. To this end we look for a non- vanishing ex- 
citonic expectation value (c^/), indicating a spontaneous 
symmetry breaking due the pairing of c electrons {f' > 0) 
with / holes {t^ < 0). Employing the normal-ordered 
representation of fermionic operators (: ... :), the Hamil- 
tonian (|T|) reads H ^Hq + Hi with 



(1) 



Ho 



with two dispersive tight-binding bands e^'^ = e"^'-^ — 
t°'-^7j^ — ij . Here e'^'-^ are the on-site energies, t'''^ are 
the nearest-neighbor particle transfer amplitudes, 7^. = 
2 cos kd for a D-dimensional hypercubic lattice, and /x 
denotes the chemical potential. Accordingly the fermionic 
operators c^^ and f^'' annihilate (create) spinless c and 
/ electrons with momentum k, respectively, and and 
nj are the corresponding particle number operators for 
Wannier site i. U parametrizes the local Hubbard at- 
traction. Note that if the c and / bands are degenerate, 
e"^ = and = t-^ , the EFKM reduces to the stan- 
dard Hubbard model [TB], whereas for = the genuine 
FKM arises [S]. In the latter case the local / electron 
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Here d-^ ~ i'^fk) pla-ys the role of the EI order parameter, 

and flkik^k,, 7 CkiCk^/ka/ki+ka-k.- Note that in Ho, the 
on-site energies are shifted by a Hartree term, 
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where {n") = jjT,k{clc^) and (nf) = ^Ek(/k/k) are 
the mean particle number densities of c and / electrons for 
a system with N lattice sites. Thus, Ho alone corresponds 
to the Hartree Hamiltonian. Vice versa, only fluctuation 
operators from the U term contribute to Hi- Evaluating 
the expectation values (...), the temperature enters into 
the calculation via the Fermi function, see [T7]- 

Following the procedure of the PRM approach [17] by 
integrating out all transitions due to T^i, the Hamiltonian 
"H can be transformed to a fully rcnormalized Hamiltonian 

k k 

+ ^(Ak:/^Ck:+H.c.), (5) 

k 

with modified parameters e^, ejj^, and Ak. Note that they 
take important correlation effects into account, which en- 
ter from the elimination procedure. 

Sclf-evidently the single particle operators have to be 
transformed in order to evaluate expectation values, i.e. 

• Ck : = ik • Ck : E ^kikks ■ "^ki /ks-^ki+ka-k -(6) 

kik2 

■Ik - = ^k'-fl- E 2/klk2,k-ki+k2 

kik2 

^ ■ CkiCk2/k-ki+k2 • ' (''') 

where ik, yk, ■ • ■ sltc also rcnormalized parameters. In the 
PRM the rcnormalization results from integrating differ- 
ence equations with initial conditions taken over from the 
original Hamiltonian ([2]) and the original single particle 

operators: e^;2 = ^k,A = ^kA^kA^ = ^' ^nd 

2/kik2k3,A(2/kik2k3,A) = 0- 

The final Hamiltonian ([S]) can be diagonalized by a Bo- 
goliubov transformation 

^ = E^k:4ck:+E^k:/k/k:+^, (8) 

k k 

where the quasiparticle energies are given by 

^c// = !k±£^ ^ '^^^'{-'-^ Wk (9) 

with Wk = m - ~^if + 4|Akni/2. 

A finite A^. signals c-/ electron coherence connected 
with a band gap that stabilizes the EI phase. Outside 
the EI phase, where Ak = 0, a band gap may also exist, 
provided that 

Ef =El-Ei=el-~el. (10) 

Therefore Ef^ < {E^^^ > 0) may be taken as indication 
that the system is in the semimetallic (semiconducting) 
regime (cf. fig. [1]). Let us emphasize that E^^ contains 
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Fig. 2: (Color online) Quasiparticle dispersions i?£ (black 
lines), (red lines) [upper panel], and (b) "order-parameter" 
functions dk [lower panel] along the major axes of the square 
lattice Brillouin zone. Results are calculated for — —1.0, 
= —0.3, and U — 2.0 at various temperatures. 

the fully rcnormalized quasiparticle energies e'^'^, not just 
the Hartree energies e'^'-^ given by eq. (|3]). 

For the 2D tight-binding band case studied in this pa- 
per, we work on a discrete set of = 24 x 24 lattice sites 
and determine all quantities with a relative error of less 
than IQ-s. 

Semimetallic region. Figure [5] (upper panel) shows 
the rcnormalized quasiparticle bands i?k and along the 
high-symmetry axes of the 2D Brillouin zone, for = — 1, 

= —0.3, and U = 2. In this case, both bands overlap 
{Eg^'^ < 0) leading to a large Fermi surface, where both 
types of quasiparticles participate. At low temperatures 
a gap opens at the Fermi surface due to the formation 
of an excitonic insulating state. Such an EI state has 
been viewed before as a BCS condensate of loosely bound 
electron-hole pairs [S]. Increasing the temperature above 
some critical temperature Tei the gap vanishes. At this 
temperature the El-semimetal transition takes place. The 
lower panel of fig. [2] displays the order parameter func- 
tion dk as a function of k. For low temperatures and k 
close to the Fermi surface, where both quasiparticle bands 
overlap, dk is strongly peaked. Otherwise dk is a rather 
smooth function of k. As a matter of course, increasing T 
above Tei, the order parameter function dk vanishes. 

In fig. [3] [panel (a)] the EI order parameter A [eq. Q] is 
shown as a function of temperature for various values of e-^ 
at C/ = 2. Clearly seen is the formation of an EI state with 
nonzero A at low temperatures. The EI state is weakened 
by lowering since the overlap of c and / electron bands 
is reduced in this case. In panel (b), the temperature- 
dependence of A is illustrated for various values of the 
Coulomb interaction U. Similar as before, the formation 
of an EI state is observed for low T. The EI region is 
decimated by lowering the c electron / hole attraction. 
One can assure oneself that the EI phase only appears 
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Fig. 3: (Color online) Temperature dependence of the EI order 
parameter A in (a) for various sX, U — 2 and in (b) for 
various U at — —2.0 where — —0.3. The solid lines give 
the variation of the gap iSg*^' (|10p (note the scale on right-hand 
ordinate. 
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in between some lower critical value Ud and some upper 
critical value Uc2 (on the semiconductor side, see below). 
In both panels of fig. [31 the solid blue lines give the varia- 
tion of the negative quasiparticle gap Eg^^ with T, where 
= —2.5 in panel (a) and U ~ 2.4 in (b). In either case, 
the variation of Ef^ for small T towards lower values goes 
along with the formation of the EI state. For U < Uci 
the bare band splitting e-^ — e'^ is somewhat reduced but 
is still negative, so we end up with a semimetallic 
situation. 

Semiconducting region. We now discuss the possible 
appearance of an EI state on the semiconductor-side of the 
schematic phase diagram shown in fig.[TJ Figure S] displays 
the temperature dependence of the quasiparticle gap Eg^'^ 
and the order parameter A for some larger values of U 
than before. The order parameter A is finite and negative 
at low T which again signals the existence of an EI phase. 
The EI phase in this region was interpreted as a BEC of 
preformed tightly bound excitons [S] . For T above the crit- 
ical temperature Tei the order parameter vanishes and no 
broken-symmetry state exists as in the semimetallic case. 
Note that E^°^ is now positive at low temperatures, which 
indicates that we have a situation with a semiconductor- 
like band structure at least up to some critical temper- 
ature Tmi, where the band gap closes. Obviously, at 
Tmi, a semiconductor-semimetal metal-insulator transi- 
tion occurs. We further note that due to the inclusion 
of correlation effects the PRM metal-insulator transition 
temperature Tmi strongly deviates from T|^j obtained by 
using the Hartree-shiftcd bare energies ej^'^ only, where 
j^{o),H _ _ j^f ^ Eq ~ Eq marks the corresponding 
Hartrce band gap. The essential question whether cxi- 
tonic bound states might possibly exist in the tempera- 
ture region Tei < T < Tmi will be investigated below. 
In this connection, in ref. [S] the authors proposed a so- 
called "halo phase" , where individual valence bond holes, 
conduction band electrons, and bound (but uncondensed) 
electron-hole pairs (excitons) should coexist. 

Excitonic resonances. — To address the possible for- 
mation of excitonic bound states above Tei, we analyze 



Fig. 4: (Color online) EI order parameter A (symbols) and en- 



ergy gaps iJg"' (solid lines), Eg'^''" (dashed lines) as functions 
of temperature for several values of the Coulomb interaction U 
where = —2.4, = —0.3. Tei and Tmi denote the criti- 
cal temperatures for the EI transition and the metal-insulator 
transition, respectively. Note that directly opposed to fig.Othe 
order parameter A decreases with increasing U in the semicon- 
ducting regime because the effective band overlap (necessary 
to establish c-f electron coherence) is now reduced on account 
of the Hartree term. 
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the frequency and momentum dependence of the dynam- 
ical excitonic susceptibility 



X(q,c.) = ((feq;6^))(c.), 



(11) 



where the symbol ((. . .)) denotes the retarded Green's 
function, and the creation operator of an electron- 
hole excitation with momentum q is defined by 6]^ = 

"^ik-Fq/k- Using the unitary invariance of expec- 
tation values, Eq. can be rewritten as 



(12) 
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Here the two-particle Green's function on the right-hand 
side is formed with H and the quantities with tilde symbols 
are the fully renormalized operators. Taking into account 
that the EI order parameter vanishes for T > Tei (Ak = 
Vk), we obtain up to order 0{U'^) 
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and T] = 0^. The coefScients F* are given by 



-lmx(q,a))/ji [arb. unit] 



-lmx(q,oi))/7t [arb. unit] 



'J^ 2;ijXlj_[_q(5;ki+qyki+q.k+q,k("-ki+q) 



~^J^ |^2a;ija;ij_l_qyk2,k:i+q,ki-k2+k+q 

kik2 

X yk2,k+q,ki-k2+k+q("-k2)("-ki-k2+k+q) 

+ 2a;ki2Sk2+qyk2,ki+q,kiyk2+q,k+q,k("-ki)("'k+q) 

+ iki+qa;k2+qyk2+q.k+q,kyki+q,k+q,k 

X Ri+q>(^k2+q) 

+ ^ki ^k2^k+q.k2+q.k22/k+q.ki+q.ki (^ki ) ("^ ) } 
x((nO-R+q))' (17) 

rkqkik2 ~^^(|Skyki.k+q,k2 P 
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Here the expectation values (nj^) = (c]{j.Cjj.}.j^ and (n^) = 
(/k/k)-?ii ^'^^ formed with the renormalized Hamiltonian 
H and can easily be evaluated due to the diagonal form 
of T-L. Note that the pole structure of the first (coher- 
ent) term of Eq. describes the continuum of particle- 
hole excitations. Of course, we have wo(0) = E^g'\ and 
in view of the form of the c and / band dispersions 
cjk(O) > wo(0) Vk 7^ 0. Therefore the possibility of q = 
excitations with positive energy indicates that the system 
is in the semiconducting regime. If one tries to determine 
the scmimetal-semiconductor boundary, i.e., Tmi, from 
the pole structure of x(q, w) this assertion is valid to lead- 
ing order only; the second and third term of (|13|) might 
lead to a shift of the lowest excitation energy in the q = 
sector. As shown below, this effect is negligible however: 
the values of Ta// derived from E'g"' are in accord with the 
results obtained from the dynamical susceptibility. 
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Fig. 5: (Color online) Imaginary part of the dynamical sus- 
ceptibility — Imx(q,a;)/7r versus frequency vertically shifted 
to show its values for different momenta along the {qx,Q)- 
direction (black solid lines). The boundaries to the particle- 
hole continuum are marked by red dashed lines. The insets 
magnify the low-frequency small-momentum region whenever 
a gap appears for the electron-hole excitations. Bare band pa- 
rameters are — —2.4, tf = —0.3. 



The imaginary part of x(q, io) reads 
-ilmx(q,a.) =i-5^FS:q<5[^-^k(q)] 



J3 X! Pkqkik2'^('^ - ^kqkika) 



kkika 



-r^qk,k2'5(c^-4qk,k2)]' (20) 

where possible non-zero excitations outside the particle- 
hole continuum point to the existence of excitonic reso- 
nances. 

Figure Odiplays our numerical results for the imaginary 
part of the exitonic susceptibility — Imx(q, w + i0)/7r as a 
function of ui for different momenta {qx , 0) between Qx ~ 
and TT (black solid lines), indicated by the scale attached 
to the ordinate axis. In all panels the red dashed line 
represents the boundaries to the particle-hole continuum 
(indicated by the strong upturns). The small bumps in 
the figures correspond to excitonic resonances, where the 
maximum most closely located to the particle-hole contin- 
uum refers to an excitonic bound state. 

In panels (a) to (d) the temperature is kept to T = 0.08. 
For U = 2.5 respectively 2.85 (which are still larger 
than Uci however), T > T^j, and the system is in the 
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semimetallic region. No excitonic resonances can be found 
for q = in this case, since Eg^^ is already negative (cf. 
fig- El)- There is a significant increase of the spectral weight 
of the finitc-q excitonic resonances by going over from 
U ~ 2.5 to U ~ 2.85. In contrast, in panel (c), where 
U — 2.9, we have Tei < T < Tmi, and the system real- 
izes a semiconductor (cf. figs. [1] and E]). Now a weak exci- 
tonic resonance is found at momentum q = which can 
more clearly be seen from the inset. For U = 3.0 > Uc2 
[panel (d)] the EI phase is not realized even for T = 
(cf. fig. HI), and again excitons can be formed for all values 
of Qx- Note that excitons with finite momentum q can 
be created in both semiconductor and semimetal cases [cf. 
panels (a) to (d)]. Their resonance positions follow from 
eqs. ^ and (fTB)) . 

In the two lowermost panels (e) and (f) the Coulomb 
interaction is fixed to U = 2.9. Increasing the temper- 
ature from T = 0.12 (e) to T = 0.16 (f) the system 
passes the semiconductor-semimetal transition. Although, 
at T = 0.12, the system is very close to the transition 
point (cf. fig. H]), excitons with zero momenta may form 
(see inset). In contrast, no q = excitons can exist for the 
temperature considered in panel (f ) . Here we observe only 
excitonic resonances with finite momenta. For the higher 
temperature case these resonances are more smeared out 
and their weight is enhanced. Therefore excitonic states 
with q 7^ can easier be occupied for this case. 

Figure [5] clearly shows that it is possible to extract 
the temperature Tmi just as well by monitoring the ap- 
pearance of excitonic resonances in the imaginary part of 
x(q, w + zO) at q = 0. 

Conclusions. — In summary, we have performed a 
detailed investigation of the two-dimensional extended 
Falicov-Kimball model by means of the projector-based 
renormalization method. Thereby we established the long- 
predicted existence of an intervening excitonic insulator 
phase at the semimetal-semiconductor transition below 
some critical Tei (see fig.H]). We derived the renormahzed 
quasiparticle band structure which shows a correlation- 
induced single-particle gap and c-f electron coherence 
in the low-temperature EI state and reflects the metal- 
insulator transition at Tmi for higher temperatures. An- 
alyzing the imaginary part of the excitonic pair suscepti- 
bility, we demonstrate that on the semiconductor side of 
this phase transition, preformed excitons with zero mo- 
mentum exist above Tei- On the other hand, excitonic 
bound states (resonances) with finite momentum may ap- 
pear on both — semiconducting and semimetallic — sides of 
the metal-insulator transition, but these excitons will not 
condense for the studied direct band gap situation. We 
therefore corroborate the scenario, suggested by Bronold 
and Fehske [5], that in the semiconducting region the EI 
phase is surrounded by an excitonic halo consisting of free 
electrons, holes and tightly bound zero-momentum exci- 
tons. Forming the EI state, the latter undergo a Bose- 
Einstein condensation state as the temperature is lowered. 



Contrariwise there is a well-defined (large) Fermi surface 
in the semimetallic regime, and the EI state can be envis- 
aged as composed of BCS-type electron-hole pairs. 
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